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Abstract: The probability density functions measured by Lewis and Swinney for turbulent 
Couette- Taylor flow, observed by Bodenschatz and co-workers in the Lagrangian measure- 
ment of particle accelerations and those obtained in the DNS by Gotoh et al. are analyzed 
in excellent agreement with the theoretical formulae derived with the multifractal anal- 
ysis, a unified self-consistent approach based on generalized entropy, i.e., the Tsallis or 
the Renyi entropy. This analysis rests on the invariance of the Navier-Stokes equation 
under a scale transformation for high Reynolds number, and on the assumption that the 
distribution of the exponent a, introduced in the scale transformation, is multifractal and 
that its distribution function is given by taking extremum of the generalized entropy with 
the appropriate constraints. It also provides analytical formula for the scaling exponents 
of the velocity structure function which explains quite well the measured quantities in 
experiments and DNS. 
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The multifractal analysis of turbulence [T]-jll| is a unified self-consistent approach for the systems 
with large deviations, which has been constructed based on the Tsallis-type distribution function 
|121I13| that provides an extremum of the extensive Reny |14| or the non-extensive Tsallis entropy 
[1211131 HB*] under appropriate constraints. The analysis rests on the scale invariance of the Navier- 
Stokes equation for high Reynolds number, and on the assumptions that the singularities due to 
the invariance distribute themselves multifractally in physical space. The multifractal analysis is 
a generalization of the log-normal model ^H1~CB|- It nas been shown jSj that the multifractal 
analysis derives the log-normal model when one starts with the Boltzmann-Gibbs entropy. 
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In this paper, we derive the formula for various probability density functions (PDF's) in fully 
developed turbulence by means of the multifractal analysis, and analyze the PDF's observed in 
three experiments. The first is the PDF of velocity fluctuations at R\ = 262 of the Taylor mi- 
croscale Reynolds number measured by Lewis and Swinney 19 in the real experiment for turbulent 
Couette- Taylor flow in a concentric cylinder system. The second is the PDF of accelerations at 
R\ = 970 obtained in the Lagrangian measurement of particle accelerations that was realized by 
Bodenschatz and co-workers [20112] by raising dramatically the spatial and temporal measurement 
resolutions with the help of the silicon strip detectors. The third is the PDF's of velocity fluctua- 
tions, of velocity derivatives and of fluid particle accelerations at R\ = 380 that was extracted by 
Gotoh et al. from the DNS of the size 1024 3 |2*2*| . 

For high Reynolds number Re 3> 1, or for the situation where effects of the kinematic viscosity 
v can be neglected compared with those of the turbulent viscosity, the Navier-Stokes equation, 
du/dt + (u ■ V)u = -V (jp/p) + vV 2 u, of an incompressible fluid is invariant under the scale 
transformation [23 El r -> Ar, u -> X a / 3 u, t -> A 1 " a/3 t and (p/p) -> A 2q/3 (p/p) where the 
exponent a is an arbitrary real quantity. The quantities p and p represent, respectively, mass 
density and pressure. The Reynolds number Re of the system is given by Re = Sui n £i n /v = 
{£ in /r]) 4/3 with the Kolmo gorov scale r\ = (j/ 3 /e) 1 / 4 [3BJ where e is the energy input rate at the input 
scale £- m . Here, we introduced Su in — \u(»+£i n )—u(»)\ with the definition of the velocity fluctuation 
Su n = \u(» + 1„) — u(»)| where u is a component of velocity field u, and l n is a distance between 
two points. The pressure (divided by the mass density) difference Sp n = \p/p(u + l n ) — p/ p(u)\ 
between two points separated by the distance £ n is another important observable quantity. We 
are measuring distance by the discrete units £ n = S n £o with S n — 2~ n (n = 0, 1,2, ■••). The 
non-negative integer n can be interpreted as the multifractal depth. 4 However, we will treat it as 
positive real number in the analysis of experiments. 

Let us consider the quantity Sx n = \x(* + £ n ) — having the scaling property \x n \ = 

\5x n /Sx \ — Si-T^ 3 . Its spatial derivative is defined by \x'\ = \iui£ n ^ Sx n /£ n cx lim^oo £ Ka /3-i 
which becomes singular for a < 3/k. The values of exponent a specify the degree of singularity. We 

see that the scale invariance provides us with Su n /6uo = S„^ 3 and 5p n /Spo = {£ n /£o) 2a ^ 3 giving, 
respectively, k = 1 for the velocity fluctuation and k = 2 for the pressure fluctuation. The velocity 
derivative and the fluid particle acceleration may be estimated, respectively, by \u'\ = limn^oo u' n 
and by |a| = limn^oo a ra where we introduced the velocity derivative u' n — Su n /£ n and the acceler- 
ation a n = Sp n /£ n corresponding to the characteristic length £ n . Note that the acceleration a of a 
fluid particle is given by the substantive time derivative of the velocity: a = du/dt + (u ■ V)u. We 
see that the velocity derivative and the fluid particle acceleration become singular for a < 3 and 
a < 1.5, respectively, i.e., \u'\ oc lim^^o ii^ 3 ^ 1 — > oo and \a\ cx lim^ n _»o ^n 2 ™^ 3 ' 1 ™> °°- 

2 Multifractal analysis 

The multifractal analysis rests on the multifractal distribution of a. The probability P^ n \a)da to 
find, at a point in physical space, a singularity labeled by an exponent in the range a ~ a + da 
is given by E]-[S] p(")(a) = [l - (a - a ) 2 /( Aa) 2 ] n/{1 ~ q) /Z^ with an appropriate partition 
function Za and (Act) 2 = 2X/[(1 — q)ln2]. The range of a is a m i n < a < a max with a m i n = 

ao — Aa, a m ax = «o + Aa. This is consistent with the relation |241 IB] P( n \a) oc 8 n ~^ a ^ that 
is a manifestation of scale invariance and reveals how densely each singularity, labeled by a, fills 

4 The multifractal depth n is related to the number of steps within the energy cascade model . At each step of the 
cascade, say at the nth step, eddies break up into two pieces producing an energy cascade with the energy-transfer 
rate e n that represents the rate of transfer of energy per unit mass from eddies of size £ n to those of size £a+i. The 
energy dissipation rate becomes singular in the limit n — > oo for a < 1, i.e., lirrin^oo e n /eo = lim n ^oo (£n/^o) a 1 —* 
oo. 
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physical space. In the present model, the multifractal spectrum f(a) is given by f(a) = 

1 + (1 — log 2 [l — (a — ao) 2 /(Aa) 2 ]. In spite of the different characteristics of the entropies, 
i.e., extensive and non-extensive, the distribution functions P^ n '(a) giving their extremum have 
the common structure. 5 

The dependence of the parameters ao, X and q on the intermittency exponent \i is determined, 
self-consistently, with the help of the three independent equations, i.e., the energy conservation: 
(e n ) — e, the definition of the intermittency exponent fi: (e 2 ) = e 2 <5~^, and the scaling relation 6 : 
1/(1 — q) = l/a_ — l/a+ with a± satisfying f(a±) = 0. The average (• • ■) is taken with P^ n \a). 

It has been shown that the probability II^ n \x n )dx n to find a physical quantity x n in the range 
%n ~ x n + dx n is given in the form 

IlW (x n )dx n = 77 s ( " ) (x n )dx n + A/7 W (x n )dx n (1) 

with the normalization dx n II^ n ' (x n ) — 1. The first term represents the contribution by the 
singular part of the quantity x n stemmed from the multifractal distribution of its singularities 
in physical space. This is given by II^ l \\x n \)dx n oc P^{a)da with the transformation of the 

variables, \x n \ = 5% a ^ 3 . Whereas the second term AIJ^ (x n )dx n represents the contribution from 
the dissipative term in the Navier-Stokes equation, and/or the one from the errors in measure- 
ments. The dissipative term has been discarded in the above investigation for the distribution of 
singularities since it violates the invariance under the scale transformation. The contribution of 
the second term provides a correction to the first one. Note that each term in Q is a multiple 
of two probability functions, i.e., the one to determine the portion of the contribution among two 
independent origins, and the other to find x n in the range x n ~ x n +dx n . Note also that the values 
of x n originated in the singularity are rather large representing intermittent large deviations, and 
that those contributing to the correction terms are small in comparison with its deviation. 

The mth moment of the variable |a; n | is given by {{\x n \ m )) = dx n \x n \ m IJ^ n \x n ) = 2 7 ^ + 
(l-2 7 £ n) ) a Km #™ where 2 7 ^ } = dx n \x n \ m AI7^(x n ), a 3 g = {2/[^Q(l + V^)]} 1/2 with 
Cg = 1 + 2q 2 (l - q)X1n2. The quantity 0-0 

Cm = a m/3 - 2Xm 2 /[9(l + C^*)] - [1 - log 2 (l + C^ a 3 )]/(1 - q) (2) 

is the scaling exponent of the mth moment of the velocity structure function {{\Su n /Suo\ m }) , i.e., 
the velocity fluctuations (k = 1). It explains the experimental results, successfully. The formula 
P| is independent of n, which is a manifestation of the scale invariance. 

We now derive the PDF, //( n )(£ n ), defined by the relation 77 (n) (£n)d£™ = n^(x n )dx n with 
the variable = x n / {(x^)) 1 / 2 scaled by the deviation {{x^)) 1 / 2 . This PDF is to be compared 
with the observed PDF's. The variable is related with a by |£„| = £„<5n Q ^ 3 ^ 2k ^ 2 with = 
[272 n) <^ C2 * + (1 - 2^ n) )a 2K }- 1/2 . It is reasonable to imagine that the origin of intermittent rare 
events is attributed to the first singular term in {Tjl. and that the contribution from the second 
term is negligible. We then have, for £* < |£„| < £™ ax , 



n (n) (tnWn = n^(x n )dx n = 



1 1-g (31n|$„/^ |) 2 
n 2n 2 X\\n8 n \ 



n/(l-q) 

(3) 



with &,,„ = UST° /3 ~ <2J2 , C ax = Un atnin/3 ~^ /2 , 77 s (n) = 3(1 - 2 7 W)/(2<n\/27rX|ln5 n |). 
On the other hand, for smaller values, the contribution to the PDF comes, mainly, from thermal 

5 Within the present formulation, the decision cannot be pronounced which of the entropies is underlying the 
system of turbulence. 

6 The scaling relation is a generalization of the one derived first in ESI W\ to the case where the multifractal 
spectrum has negative values. 
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fluctuations or measurement error. It may be described by a Gaussian function, i.e., for |£„| < £*, 



fl (n) (£„)<£„ 



n& n) [x n ) + Aii w (xn 



(n)/ 



dx„ = 77 s (n) exp | - - 



[(Cn/C) 2 -!]}^. (4) 

This specific form of the Gaussian function is determined by the condition that the two PDF's J3j) 
and (QJ should have the same value and the same slope at £* which is defined by £* = In^n" ^~ < ' 2k ^ 2 
with a* being the smaller solution of C2«/2 — Ka/3 + 1 — f(a) = 0. It is the point at which fiW (£*) 
has the least n-dependence for large n. 

With the help of the second equality in we obtain AII^ (x n ), and have the formula to 

evaluate 7 i n) in the form 2 7 ^ } = (k^ ] - L& } ) / (l + - L n) ) where 



(«) - 



LW = 



3 <5. 



K(m+l)Q*/3-C2 K /2 



Ky/27tX| lni5„ Jo 



z m exp < — 



1 



i + -/V) 

K 



,, c Kmoc* /3 

3 o„ 



Ky27rX| In <5„ | 

C If* — cK(a««-tt")/3 
Smm/ — 



1 - 



1 - q (31n|z/zo|)' 



z 2 -l) 

!»/(!-?) 



2K 2 X|ln<5„| 



(5) 
(0) 



?o = Cn,o/C = Sn {a °~ a " )/3 . Now, the PDF for the variable 



with z* liri 

scaled by its own deviation, given by and Q, is completely determined by the intermittency 
exponent fi and the multifractal depth n which gives a length scale l n . The PDF's for velocity 
fluctuations and for derivatives are given by © and l@J with k = 1, whereas those for pressure 
fluctuations and for fluid particle accelerations with k = 2. The PDF for energy dissipation rates 
is given with k = 3. 



3 Analysis of experiments 

The PDF's for velocity fluctuations measured by Lewis and Swinney ^Hj at R\ = 262 for turbulent 
Couette- Taylor flow and those extracted by Gotoh et al. from his DNS data [22] at R\ = 380 are 
shown, respectively, in Fig. ^(i) and in Fig.^-(ii). They are analyzed with the derived formulae 
© and Q with k = 1 HHE|. 

For Fig.^(i), we adopted the reported value fi = 0.28 ^5] to calculate the parameters q = 0.471, 
«o = 1.162 and X = 0.334. The dependence of n on r/rj is extracted as @ 

n = -1.019 x log 2 r/rj + 0.901 x log 2 Re (7) 

with Re = 540 000. The Reynolds number is estimated with t m « 119.32 cm and the Kolmogorov 
scale r\ ss 0.006 cm. The energy input scale i- m is given by the size of experimental apparatus 
2ir x 19.00 cm The definition of the number of steps n within the energy cascade model is 
given by n — — log 2 (r/^i n ) for the eddies whose diameter is equal to r. By putting r = £ n , this 
gives us the relation between n and n in the form 

n = n-log 2 {e /£ ia ). (8) 

For Fig. n-(ii), we extracted the value fi — 0.240 by analyzing the measured scaling exponents 
( m of velocity structure function with the formula J2J, which gives the values q = 0.391, ao = 1-138 
and X = 0.285. Through the analyses of the PDF's for velocity fluctuations, we extracted the 
formula for the dependence of n on r/rj: [71 |H] 



n = -1.050 x log 2 r/?7+ 16.74 (for £ c < r), 
n = -2.540 x log 2 r/rj + 25.08 (for r < l e ). 



(9) 
(10) 
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This shows that the inertial range is divided into two scaling regions separated by the characteristic 
length £ c /rj = 48.26 which is close to the Taylor microscale A/77 = 38.33 of the system. The equation 
© is consistent with the picture of the energy cascade model in which each eddy breaks up into 
two pieces at every cascade steps, whereas ijTHt indicates that, for r < £ c , each eddy breaks up, 
effectively, into 1.33 [Hj pieces at every cascade steps. This fact may be attributed to a manifestation 
of structural difference of eddies, which can be checked by visualizing DNS eddies. Actually, one 
observes that DNS eddies with larger diameters than Taylor microscale A have rather round shapes, 
whereas eddies with smaller diameters compared with A have rather stretched shapes |28| . The 
energy input scale for this DNS is estimated as the longest scale available in the lattice with cyclic 
boundary condition, i.e., l- m jr] = n/r] ~ 1220 with 77 w 0.258 x 10 -2 [22]. With this value of £- 1TL , 
© with © and JEJ gives the number of steps n within the energy cascade model for the DNS. 

The analysis of the PDF for velocity derivatives reported by Gotoh et al. [221 are performed 
with /j = 0.240. We chose the value n = 23.1 (n = 17.4). The corresponding length r = l n 
is calculated by 1|10JI to give rjr\ = 1.716. This length may give us an estimate for the effective 
shortest length in processing the DNS data to extract velocity derivatives. Note that it is about 
the same order of the mesh size Ar/77 = 2tt/(1024 x rf) w 2.38 |22] of the DNS lattice. 



(i) 




5 10 10 20 30 



Figure 1: PDF's for longitudinal velocity fluctuations (closed circles; symmetrized by taking av- 
erages of the left and the right hand sides data) observed by (i) Lewis and Swinney at R\ = 262 
(Re = 540 000), and by (ii) Gotoh et al. at R\ — 381. For the experimental data, the distances 
r/ri = £ n /T) are, from top to bottom: (i) 11.6, 23.1, 46.2, 92.5 208, 399, 830, 1440; (ii) 2.38, 4.76, 
9.52, 19.0, 38.1, 76.2, 152, 305, 609, 1220. For the theoretical PDF's (solid line), from top to 
bottom: (i) fi = 0.28, n (n) = 14 (10.7), 13 (9.7), 11 (8.7), 10 (7.7), 9.0 (6.6), 8.0 (5.6), 7.5 (4.6), 
7.0 (3.8); (ii) fi = 0.240, n (n) = 21.5 (15.8), 20.0 (14.3), 16.8 (11.1), 14.0 (8.31), 11.8 (6.11), 10.1 
(4.41), 9.30 (3.61), 8.10 (2.41), 7.00 (1.31), 6.00 (0.31). For better visibility, each PDF is shifted 
by — 1 unit along the vertical axis. 

The PDF's for fluid particle accelerations measured by Bodenschatz et al. at R\ = 970 |20II21| 
and those extracted out from the DNS data by Gotoh et al. at R\ — 380 22J are shown, respectively, 
in Fig. EK 1 ) and in Fig. |2J-(ii), on log and linear scale [101 111) . They are analyzed with the derived 
formulae © and J3} with k = 2 jEHH- 

For Fig.[21-(i), we determined the value n = 17.1 for this experiment by substituting the reported 
value 7.1 cm for i$ and the spatial resolution 0.5 /xm of the measurement for £ n into its definition, 
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n = \og 2 (£o/£ n ). The intermittency exponent \l = 0.250 is extracted by the analysis of the 
experimental PDF with the derived theoretical formula |101 lllj . Then, we have the values of 
parameters: q = 0.413, «o = 1-144 and X = 0.297. The flatness of the PDF turns out to be 
= ((a^))/((a^)) 2 = = 61.3 which is compatible with the value of the flatness ~ 60 

reported in |20l I21| . 

For Fig. 0-(ii), with \i = 0.240 for this DNS, we have the values of parameters: q — 0.391, 
a = 1.138, X = 0.285. The analysis of the PDF obtained by the DNS gives the value n = 18.3 
(n = 12.6). Substitution of this value into (fTTTf) gives the corresponding characteristic length 
r/rj = 6.36 This may be the effective minimum resolution in cooking the DNS data to distill 
accelerations. 




Figure 2: PDF's for fluid particle accelerations measured by (i) Bodenschatz et al. at R\ = 970 
(open circle), and by (ii) Gotoh et al. at R\ = 380 (closed circle), plotted on (a) log and (b) linear 
scale. The experimental data points both on the left- and right-hand sides are plotted altogether. 
For the theoretical PDF's (solid line), (i) fj, = 0.250 and n = 17.1, and (ii) fx = 0.240 and n = 18.3 
(n = 13.6) (solid line). 



4 Discussions 

It was shown that the various experimental PDF's in turbulence are analyzed successfully with the 
formulae J2J and (@J derived by the multifractal analysis. 

From the above analyses, we see that the contribution of thermal fluctuation and/or measure- 
ment error to PDF's is restricted to smaller values, i.e., < £*. In the case of the PDF's for ve- 
locity fluctuations, ^ = 1.10 - 1.32 (a* = 1.08) for Lweis and Swinney QH, and = 1.01 - 1.39 
(a* = 1.07) for Gotoh et al. $2\. As for the PDF for velocity derivatives by Gotoh et al. [22, 
£* = 0.982 (a* = 1.07). In the case of the PDF's for fluid particle accelerations, £* = 0.565 
(a* = 1.01) for Bodenschatz et al. HUE], and £* = 0.551 (a* = 1.005) for Gotoh et al. 
Within the present approach, the intermittent large deviations £„ > £* are a manifestation of the 
multifractal distribution of singularities a due to the scale invariance of the Navier-Stokes equation 
for Re > 1. 

There is no room to incorporate into the present multifractal analysis the energy input scale 
£i n and the "system size" £q. The former is necessary to determine the number of steps n in 
the energy cascade model. Once t- m is determined by investigating the structure of experimental 
apparatuses, the relation between n and the multifractal step n is given by (jHJ. Since main part of 
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the multifractal analysis rests on the scale invariance, the size of the system under consideration 
is assumed to be infinite, and therefore, the length £q may not have important physical meaning. 
Actually, the empirical equation J7J) extracted from the experimental PDF's for velocity fluctuations 
by Lewis and Swinney gives £$ ~ 877 cm which is large compared with the largest size l- m w 
119.32 cm of the experimental apparatus. For Gotoh's DNS, the empirical equation © gives 
@o/v ~ 63000 which is larger than the largest size £i n /r] ~ 1220 of the DNS lattice. On the other 
hand, in the analysis of Bodenschatz's experiment, the assignment of £q to the integral length scale 
7.1 cm gives us reasonable value n. It is worthwhile to note here that PDF's derived within the 
multifractal analysis seem to be sensitive to the characteristic lengths such as the distance of two 
measuring points, the space resolution in measurement and the mesh size of DNS. How to put 
the information of characteristic lengths of experimental apparatus into the multifractal analysis 
is one of the important future problems. It may be resolved when one succeeds to reveal the 
dynamical foundation underlying the basis of the multifractal analysis starting an investigation by 
the Navier-Stokes equation with the energy input term. 

The authors would like to thank Prof. R.H. Kraichnan and Prof. C. Tsallis for their fruitful 
and enlightening comments with encouragement, and are grateful to Prof. E. Bodcnschatz and 
Prof. T. Gotoh for the kindness to show their data prior to publication. 
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